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Проведен анализ динамических характеристик процесса механической обработки. Построена математиче- 
ская модель привода главного движения пилорамы. 
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Введение. Важнейшим звеном в лесопильном производстве являются лесопильные рамы. Их су- 
щественным недостатком являются неуравновешенные силы инерции подвижных масс кривошип- 
но-шатунного механизма, вызывающие вибрацию. Для анализа явлений, происходящих в меха- 
низме пилорамы и оценке технологического процесса, представляется возможным использование 
теории регулярных колебательных систем [1], основываясь на теории цепочек, позволяющей 
аналитически исследовать динамические характеристики систем с большим числом степеней сво- 
боды, базируясь на анализе одного структурного элемента [2-6]. 

Исследование процесса силового взаимодействия при распилке лесоматериала на пило- 
рамах неизбежно сводится к построению математической модели механизмов, формообразующее 
движение которых и приводит к его образованию. Одна из парциальных систем пилорамы совер- 
шает неравномерное программное движение, имеется в виду кривошипно-шатунный механизм 
привода рамы. Кроме того, на этот привод накладываются нежелательные колебания, вызванные 
кинематическим возбуждением. Согласно многочисленным работам по этой теме, особенно зна- 
чительные динамические ошибки возникают из-за виброускорений [1, 6]. Одной из основных за- 
дач в снижении виброактивности, а соответственно и уровня акустической эмиссии исследуемого 
процесса, является обеспечение требуемого закона движения (гармонического) инструмента. Та- 
ким образом, особое значение представляет исследование формообразующих движений. При ис- 
следовании динамики упругих элементов машин и механизмов, как правило, предполагается, что 
влияние вынуждающих колебания устройств одностороннее, т.е. отсутствует обратное влияние 
упругой подсистемы на источник энергии. 

Целью настоящей работы является создание математической модели привода главного 

движения пилорамы, которая и позволит решить все указанные выше задачи. 
Построение математической модели привода главного движения пилорамы. Рассмот- 
рим механизм привода пилорамы, представляющий собой последовательно соединённую ремён- 
ную передачу, и кривошипно-шатунный механизм, который реализует кинематическую функцию 
положения ведомого звена. 

Рассматриваемый механизм относится к классу устройств, преобразующих вращательное 
движение приводного вала в неравномерное движение рабочего органа. В таком случае зависи- 
мость, связывающая положение ведомого звена с углом поворота приводного электродвигателя, 
является нелинейной. Анализ кинематической схемы исследуемого механизма позволяет сделать 
вывод, что динамическая модель пилорамы может быть классифицирована как неголономная сис- 
тема из-за наличия ремённой передачи. Однако при этом нужно учитывать возможность проявле- 
ния параметрических резонансов, так как приведённый момент инерции механизма также являет- 


ся функцией ф. 
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В этом случае динамическая модель исследуемого объекта представлена в виде расчетной 
схемы (рис.1). 





Рис.1. Расчётная схема привода пилорамы 


Так как звенья этого механизма совершают сложное плоское движение в вертикальной 
плоскости, необходимо учитывать работу сил тяжести. Кроме того, в механизме имеется звено, 
совершающее плоскопараллельное движение. Для таких звеньев, согласно [4], инерционные ха- 
рактеристики могут быть приведены к звеньям, совершающим вращательное и поступательное 
движения, т.е. к кривошипу и раме. Используя метод замещения масс, получим: 

В р 
В И Ат=т,„.— 

Электромеханическая система не может быть описана адекватной моделью без учёта 
инерционности процессов, протекающих в двигателе [3]. Уравнения, описывающие движения 
асинхронного электродвигателя, имеют вид: 
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где 2, — число пар полюсов электродвигателя; Ч’, — проекция потокосцепления обмотки стато- 
ра на ось х; Ч, - проекция потокосцепления обмотки статора на ось у; Ф,, - проекция 


потокосцепления обмотки ротора на ось х; Ч, — проекция потокосцепления обмотки ро- 


2у 
тора на ось у; Л -— сумма момента инерции ротора электродвигателя и шкива; <, - угло- 
вая частота вращения поля статора; © - угловая частота вращения ротора; ‹«, - угловая 


частота. 
Таким образом, исследуемая электромеханическая система включает в себя электродвига- 
тель, описываемый системой нелинейных дифференциальных уравнений первого порядка, а так- 
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же системой уравнений, описывающих механическую часть привода. Для подобных задач обычно 
используют особую форму уравнений Лагранжа П рода с «лишними» координатами (их ещё на- 
зывают уравнениями Феррерса) [1, 2, 6]. Для рассматриваемого случая имеем: 


ня ОТ, ЭР _ ‚+ а, (Е=1....,5); 


@ 04, 04, “аа, > 
Уа,4, +а, =0, (Еаыя) 
К=1 


Обобщённые координаты выбираем, пользуясь рекомендациями, приведенными в работе 
[6]. Задаёмся абсолютной координатой, соответствующей перемещению в начале кинематической 
цепи ф, =4,, И далее вводим координаты, двигаясь по кинематической цепи. Выпишем основные 
кинематические соотношения: 





ф=О-1=9,; 
ф. =фи +4, =ай +4»; 

ф; =$, 1+9, +9; =Чй +4, +4; 

У=П(ф,) = п (ай +4, +4; т 

Е =7+4.; 2=4,+4., 
где © - угловая частота вращения приводного электродвигателя; ф,, ф,, ф, - угловые коорди- 
наты соответствующих сечений в абсолютном движении; У — абсолютная координата мас- 
сы т. Последнее замечание означает, что в качестве обобщённых координат, за исключе- 
нием а,, приняты относительные координаты, отвечающие за деформацию упругих элемен- 

тов. 

Уравнение связи запишем исходя из условия, что Х является функцией положения [1, 6]. 

П(ад 4 +4:)-2=П($,)-2=0. 

Продифференцируем это выражение по времени: 

Пай +049, +19,-й=0. 

В качестве «лишней» координаты примем У =д, . Чис- 
ло степеней свободы исследуемой системы Н=4, число и 
«лишних» координат п =1. 

В связи с тем, что рама с пилами движется в верти- 
кальной плоскости (рис.2), её перемещение удобно предста- 
вить в следующем виде: 5 =Ё+Р. 

5=Ё.с056—А-с05ф. 5 б | 
Общий катет прямоугольников ОАС и АВС определяет- 








} . [д К 
ся равенством вида: зтб = вет) ‚ где с= г. 


ЕР у 
п($)=5=2-0085+ -зтф- РВ; м 


п($) =. Е. - 05? ф+А-5шф- УД = Рис.2. Кривошипно-шатунный 


механизм пилорамы 
После преобразования имеем: 


а ааа 
аф аф 
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зе ее ааа 
аф\ аф аф 
2 2 
а - ыы -5Шб=-с.с0$ф; 
аф аф 
425 1 [, (31$) 
а 3 Е) : =} 
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п"(6)=2 Оо О соз5 ее — В созф; 
аф аф аф 
. . 2 . 2 
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Проведём оценку спектра координаты У с целью выявления частот возбуждения меха- 
низма механической частью привода. С этой целью проведем расчёт с использованием возможно- 
стей системы Масаа (рис.З3): 


В,= 0.3 - радиус кривошипа; 
Г. := 2 - длина шатуна; 


п := 325; 
© := 2 я = 34.034 - угловая скорость кривошипа; 
В 
= т = 0.15 - отношение длин кривошипа и шатуна; 


к:= 1. 1000; Щ := 0.0005; Фк = О; дк = вв _ф 9). 


пк := (| Е С2соз( к) + В-з(% ‚) = 


51 ($ к) 


Шк := ый + В-соз(% к)! 


: ь 2 | 2 
И |) м. ее) 
а соз(5 соз(5к) ь ь соз(5 } г 
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Рис.3. Соотношение выходного смещения нелинейного звена, 
его скорости и ускорения от угла поворота ведущего звена 


Проведём Фурье-анализ с использованием системы Ма{саа (рис.4). 
Как видно из анализа полученных графиков, наличие кратных гармоник просматривается, 


в наличии фазовые сдвиги, отличные от м . Последнее означает, что можно ввести некоторые 


упрощения в систему уравнений динамики механизма. 


Ф 
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Рис.4. Фурье-спектр выходной координаты нелинейного элемента 


Для составления математической модели механической части привода пилорамы выража- 
ем кинетическую и потенциальную энергию через принятые обобщённые координаты, включая и 


«лишние»: 
27 =? +.Л,ф2 +162 +т- 1; 
2Р=С, (ф, -ф,) С, (ф. -$,) +С, .(2-У). 
Диссипативная функция системы (функция Релея): 
: . \2 . . р) : 
2Ф=А ($, -$,) +^, ($, -$,) +17”. 


Выражение для виртуальной работы может быть записано в виде: 
бА = Ма-54.-Е„.-67-Ь, 4,64, -В.2 87 , 


где Ь, — приведённый коэффициент линейного демпфирования привода и ведомого звена, полу- 
ченный путём эквивалентной линеаризации диссипативных сил [1, 4, 6]; Б, — коэффициент 


вязкого трения при резании. 
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Запишем выражение кинетической энергии с учётом кинематических соотношений: 
2Т = Л4? +7, (41 +4.) +Л (ав +4, +4.) +т-(4. +4.) ; 
27 =Л 4: +7, (и +2494, +4: )+/ (ий +214, -(4,+4)+4 +4 +24.4,)+т. (4: +244; +4: ); 
27=4:(Л+лй чл +4» (Л, ++ Л, +244, (Ла + Ла) +244, + 
+24, + та: +2т4 4; + тах. 
Исходя из матричного представления кинетической энергии системы 21 а” М. а И ТОГО, 


что коэффициенты при обобщённых скоростях постоянны, инерционные коэффициенты прирав- 
ниваем к соответствующим коэффициентам квадратичной формы, представляющей кинетическую 
энергию системы, получим матрицу инерционных коэффициентов в виде: 


Л+Ле+ЛЕ Ла+ла Ла 0 


а ай МР, > НО 
АЙ ы 3» 30 
0 0 О т 


Согласно работам [1, 4] необходимо рассматривать отсчёт потенциальной энергии от по- 
ложения устойчивого равновесия Р(О,0,..,0)=Р|) =0. Кроме того, считаем, что все связи ста- 
ционарны, процесс резания рассматривается не как составляющая механизма, а как внешнее воз- 
действие. Если же принять, что передаточное отношение не зависит от скорости звеньев, матри- 


ца жёсткостей аппроксимирующей системы дифференциальных уравнений исследуемого объекта 
может быть принята постоянной [4]. 


2 2 2 
2Р=С($,.-$,) +С, ($, -$,) +С,(2-У) ; 
. 2 : р 2 2 
2Р=С (4, (1 -1)+4,) +С, (Чи 4 + 43 —414>) + С, (4, + 44 -а;) ; 
: 2 : 
Ре С, -(& -1) + 4,9>С, (1 -1)+9>-С, а: Са, С 
Исходя из матричного представления потенциальной энергии системы 2Р=а’-С-а ито- 


го, что жёсткости связей постоянны, коэффициенты матрицы жёсткостей получим приравнивани- 
ем соответствующих коэффициентов квадратичной формы, представляющей потенциальную 
энергию системы, в виде: 


са < 00 
ое © оо 
0 0 ©. 0 

0 0 о С, 


В выражение потенциальной энергии вошла координата Х’, но хотя она посредством не- 
линейной функции п($,) жёстко связанна с координатой а, , изменение потенциальной энергии 


исследуемой системы определяется перемещением инерционных масс в вертикальной плоскости 
[1, 4, 6]. 
Представим функцию демпфирования в матричном виде: 


2Ф= 1 (6, -$,) +1, ($,-5,) +1, (2-7) 


2Ф- 1 (4, (1-1) +4.) +1, (44 +4, +4, -44, + (4,+4.-4:).; 


2 
, 
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2Ф=й-й (1-1) +42 1+4 № +1, -4; +244, В (1-1). 


Исходя из представления матрицы демпфирования в виде 2Ф=а'.Н.а, получим матри- 


цун: 
в (4-1) #в.(- 0 0 
и_| №0 й, оо 
0 о ло 
0 0 ол, 


Обозначим диссипативные силы К, ‚ сохранив нумерацию, принятую для упругих элемен- 
тов. В данном случае А, — момент диссипативных сил в ремённой передаче, А, и А, - момент 
диссипативных сил, связанных с внутренним трением в валах. Моментами сил трения в опорах 
пренебрегаем [1, 6]. Дальнейшие действия сводятся к записи системы дифференциальных урав- 
нений движения и исключению множителей Лагранжа [1, 6]. Полное число уравнений равно 
Н+п=5. В первой части, помимо обобщённых сил, стоит слагаемое Ла,,, так как п=1. Запи- 


шем уравнения, устанавливающие связь «лишних» и независимых координат. 
Определим коэффициенты уравнений дополнительных связей: 


У=а. =П(5)=П(а4 +4, +9: ); 
4; =П ($, ) 44 +1 ($,)- 4, +П ($,)- 45 или п ($,)- 44 +П ($, ) 45 +П ($,)-4: 4: =0. 


Общий вид уравнения связи: 
Н+п 


У аа, +а1=0, 
К=1 
где Н +п=5. 
Сопоставим полученный результат с уравнением связи в общем виде [1,6]. Учитывая, что 
индекс коэффициента соответствует обобщённой скорости, запишем: 
а,4, + а,4> + а34; +а.4. +а;4; +а1=0. 
Из последнего равенства следует: 
а =П ($,)-:1; а, =П($,); а, =П ($,); а.=0; а, =-1. 
Определение обобщённых сил проводим, составляя выражение суммы работ на виртуаль- 
ных перемещениях: 


бА = Ма - 94, - К, - 64, - В, 64; —та-5(4. +9; )— Рте2`6(44 +4; }; 
О =Ма; О, =—№,; О, =-В,; О, =-та — Еге2; О; =-те — Етер. 
Окончательно систему уравнений вынужденного движения механизма можно записать в 


матричном виде, выделив пятое уравнение, которое предназначено для определения множителя 
Лагранжа Л: 


Ма+ Ла, Ма+АП\ 
аа ЕЛ —А, АП’ 
Ма+На+ Са= 2 — м 

—К, + Ла, —К +АП 
—та — Еге2 + Ла. —тя — Ете2 


та, +т4; = -тя — Егер -— Л. 


Последнее уравнение приведённой системы используем для определения множителя Ла- 
гранжа: 





АЛ =-тд, — тд; — те — Егез. 
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Окончательно систему уравнений привода пилорамы получаем в виде матричного равен- 


ства 
(М.р’+Н.р+С)-а=Е(Х,1); 
Ма-+Ап\ | | Ма-п\-(т4, + тд; + тв + Еге?) 
аль -® + АП’ |_| -В, -П'.(т4. + тё; + тв + Еге? ) 
ы —Ю, +ЛАП’ —К, -П’. (т4, + тф; + та + Ете2) 
—те — Ете? —те — Еге? 
Ма —П\ - (та + Еге2) й (4. +4.) 1 - (44 +45) 
:: ка т. 414% нА 4+4 
—Ю, -П - (те + Ее?) 94 +9; 94 + 4$ 
—те — Еге2 0 0 
Проведём преобразования вектора, определяемого «лишней» координатой: 
4, =П”. (1-4) +П’. пт; 
.. .. л 2 
1 - (4. +4.) (1у.а) 
р ни 
+ . 
И т = тп’. [у 4, + ТП’. п -а+ П.П" И : 
94 +45 (1.4). 
й 0 
где 
р 1то 
1 и: 80 
= Е ш=|. о 
1 |110 
0 0000 


Матрица инерционных коэффициентов, очевидно, будет иметь вид: 
М=М+[.т.П'+ а.т.0’.П". 
Тогда система уравнений, описывающая механическую часть привода, может быть пред- 
ставлена матричным уравнением вида: 


(м4) 
т ь5 2 
ма+На+Са=Е (1) -тП’.П” (т-4) = Езат. 
(1.9) 
0 


Таким образом, нами получено матричное нелинейное дифференциальное уравнение вто- 
рого порядка относительно выбранных обобщённых координат. Для дальнейшего исследования 
желательно иметь систему дифференциальных уравнений первого порядка. С этой целью прове- 
дём следующие преобразования: матрицы М,Н,С рассматриваем как составляющие блочных 

Н М С Мим Езат 


матриц: 
ЧЕ чт 
м мм Мм -мМ МУ 


где ММ - нулевая матрица 4х4; МшУ - нулевой вектор 4х1. 


7 
7 


> 
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Система дифференциальных уравнений первого порядка может быть записана в виде: 
В.Х+К.-Х=Е: О, 
где Е - единичная матрица 4х4. 
Проводим проверку сделанных преобразований: 


М Ма С м с 
м мм а| |№м ма |№У 
В результате получаем: 
Н.а+а:М+С-4а=Езат; 
М.а-М-а= ММ. 
Этот результат означает, что нами получена система уравнений, описывающая поведение 
станка в переменных состояния. Собственная матрица объекта получается в виде: 








А=-В".К, 
матрица управления в виде: 
В=В".Е, 
где 
п ($) --2.9 6+ В. совф-Е-с-50® 305+ В-с086; 
аф с0$9 
. . 2 . У 
о с $ —с:с05ф |- [с05б. оо — В. с0$ф. 
соб с0$б с0$6 


Заключение. Полученная система уравнений привода в дальнейшем может быть использована 
для анализа соотношения спектра колебаний, возбуждаемых приводом, как за счёт его динамиче- 
ских свойств, так и за счёт существенной нелинейности закона движения рабочего органа и соб- 
ственных частот пилы. Кроме того, она может быть использована и при анализе вынужденных 
движений пилы. Таким образом, поставленная цель настоящей работы реализована, создана ма- 
тематическая модель привода главного движения пилорамы. 
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Тре аупате спагасейеНс$ о! пе тесвагка! орегайоп аге апгиугеа. Те таетайса! тодЕ! о! {пе атмуе о! {те 
рои’ег-заи/ Бепсв рипсра/ тоуетепе 15 БийЕ. 
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